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Abstract 

We prove a global in time existence theorem for the initial value prob- 
lem for the Einstein-Boltzmann system, with positive cosmological con- 
stant and arbitrarily large initial data, in the spatially homogeneous case. 



' in a Robertson- Walker space-time. 

o ■ 

1 Introduction 



Jh^' In the mathematical study of General Relativity, one of the main problems is 

to establish the existence and to give the properties of global solutions of the 
Einstein equations coupled to various field equations. The knowledge of the 
global dynamics of the relativistic kinetic matter is based on such results. In 
■ the case of CoUisionless matter, the phenomena are governed by the Einstein- 

I Vlasov system in the pure gravitational case, and by this system coupled to other 

fields equations, if other fields than the gravitational field are involved. In the 
collisionless case, several authors proved global results, see P, ^Hl for reviews, 
|17|. [21] and [S] for scalar matter fields, also see [20], [21] for the Einstein- Vlasov 
system with a cosmological constant. Now in the case of Collisional matter, 
the Einstein- Vlasov system is replaced by the Einstein-Boltzmann system, that 
seems to be the best approximation available and that describes the case of 
instantaneous, binary and elastic collisions. In contrast with the abundance of 
works in the collisionless case, the literature is very poor in the collisional case. 
If, due to its importance in collisional kinetic theory, several authors studied 
and proved global results for the single Boltzmann equation, see [5] , [^ , for 
the non-relativistic case, and [3|, for the full relativistic case, very few authors 
studied the Einstein-Boltzmann system, see [H] for a local existence theorem. 
It then seems interesting for us, to extend to the collisional case, some global 
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results obtained in the collisionless case. This was certainly the objective of the 
author in ^21 and ^1], in which he studied the existence of global solutions 
of the Einstein-Boltzmann system. Unfortunately, several points of the work 
are far from clear; such as, the use of a formulation which is valid only for the 
non-relativistic Boltzmann equation, or, concerning the Einstein equations, 
to abandon the evolution equations which are really relevant, and to concentrate 
only on the constraint equations, which, in the homogeneous case studied, reduce 
as we will see to a question of choice for the initial data. 

In this paper, we study the collisional evolution of a kind of uncharged 
massive particles, under the only influence of their own gravitational field, 
which is a function of the position of the particles. 

The phenomenon is governed, as we said above, by the coupled Einstein- 
Boltzmann system we now introduce. The Einstein equations are the basic 
equations of the General Relativity. These equations express the fact that, the 
gravitational field is generated by the matter contents, acting as its sources. 
The gravitational field is represented, by a second order symmetric 2-tensor of 
Lorentzian type, called the metric tensor, we denote by g, whose components 
5(q,/3, sometimes called "gravitational potentials", are subject to the Einstein 
equations, with sources represented by a second order symmetric 2-tensor we 
denote T^p, that summarizes all the matter contents and which is called the 
stress-matter tensor. Let us observe that, solving the Einstein equations is 
determining both the gravitational field and its sources. In our case, the only 
matter contents are the massive particles statistically described in terms of their 
distribution function, denoted f, and which is a non-negative real valued 
function, of both the position and the momentum of the particles, and that 
generates the gravitational field g, through the stress-matter tensor Taf3- The 
scalar function f is physically interpreted as " probability of the presence density" 
of the particles during their collisional evolution, and is subject to the Boltzmann 
equation, defined by a non linear operator Q called the "Collision Operator". 
In the binary and elastic scheme due to Lichnerowicz and Chernikov(1940), we 
adopt, at a given position, only 2 particles collide, in an instantaneous shock, 
without destroying each other, the collision affecting only their momenta which 
are not the same, before, and after the shock, only the sum of the 2 momenta 
being preserved. 

We then study the coupled Einstein-Boltzmann system in (g,f). The system 
is coupled in the sense that f, which is subject to the Boltzmann equation 
generates the sourcesTc^ of the Einstein Equations, whereas the metric g, which 
is subject to the Einstein equations, is in both the Collision operator, which is 
the r.h.s of the Boltzmann equation, and in the Lie derivative of f with respect 
to the vectors field tangent to the trajectories of the particles, and which is the 
l.h.s of the Boltzmann equation. 

We now specify the geometric frame, i.e. the kind of space-time we are look- 
ing for. An important part of general relativity is the Cosmology, which is the 
study of the Universe on a large scale. A. Einstein and W. de Sitter introduced 
the cosmological models in 1917; A. Friedman and G. Lemaitre introduced the 
concept of expanding Universe in 1920. Let us point out the fact that the 
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Einstein equations are overdetermined, and physically meaning symmetry 
assumptions reduce the number of unknowns. A usual assumption is that the 
spatial geometry has constant curvature which is positive, zero or negative, re- 
spectively. Robertson and Walker showed in 1944 that "exact spherical symme- 
try about every point would imply that the universe is spatially homogeneous" , 
see 20], p. 135. We look for a spatially homogeneous Friedman-Lemaitre- 
Robertson- Walker space-time, we will call a " Robertson- Walker space-time", 
which is, in Cosmology, the basic model for the study of the expanding Uni- 
verse. The metric tensor g has only one unknown component we denote a, which 
is a strictly positive function called Cosmological expansion factor; the spatial 
homogeneity means that a depends only on the time t and the distribution func- 
tion f depends only on the time t and the 4-momentum p of the particles. The 
study of the Einstein-Boltzmann system then turns out to the determination of 
the couple of scalar function (a, f). 

In the present work, we consider the Einstein Equation with cosmological 
constant A. Our motivation is of a physical point of view. Recent measurements 
show that the case A > is physically very interesting in the sense that one 
can prove, as we will see, that the expansion of the Universe is accelerating; 
in mathematical terms, this means that the mean curvature of the space-time 
tends to a constant at late times. For more details on the cosmological constant, 
see [131. 

Let us now sketch the strategy we adopted to prove the global in time exis- 
tence of a solution (a, f ) of the initial values problem for the Einstein-Boltzmann 
system for arbitrarily large initial data (ao,/o) at t = 0. In the homogeneous 
case we consider, the Einstein Equations are a system of 2 non-linear o.d.e for 
the cosmological expansion factor a; the Boltzmann equation is a non-linear first 
order p.d.e for the distribution function /. 

In a first step, we suppose a given, with the only assumption to be bounded 
away from zero, and we give, following Glassey, R.T., [?], the correct formulation 
of the relativistic Boltzmann equation in f, on a Robertson- Walker space-time. 
We then prove that on any bounded interval / — [ioi^o + T] with S 
T € , the Cauchy problem for the Boltzmann equation has a unique solution 
/ € C[I; L^iR^)], where ^2(1^^) is a weighted subspace of L^(R^), whose weight 
is imposed by the expression of the sources terms Tap of the Einstein equations. 
We follow the method developed in ^H] that prove a global existence of the 
solution / e C[[0, +oo[; for the Cauchy problem for the Boltzmann 

equation , but here, the norm of i2(I^'^) could allow us to prove the existence 
theorem only on bounded intervals [to, to + T], and this was enough for the 
coupling with the Einstein Equations. 

In a second step, we suppose f given in C[I; L2{M-'^)], and we consider the 
Einstein equations in a, that split into the constraints equations and the evo- 
lution equation. The constraint equations contain the momentum constraint 
which is automatically satisfied in the homogeneous case we consider, and the 
Hamiltonian constraint that reduces to a question of choice for the initial data. 
The main problem is then to solve the evolution equation, which is a non-linear 
second order o.d.e in a. We set: e = i, 6* = 3^ (where a = ^), and we 
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prove that, the Einstein evolution equation in a is equivalent to the non-linear 
first order system in (e, 9) defined by e = — ^ and the non-linear first order 
Raychaudhuri equation in ^, {6 is called the "Hubble variable"), and for which 
we deduce that there could exist no global solution in the case A < 0. In the 
case A > 0, we deduce by applying standard theorems to the first order system 
quoted above, the existence of the solution a of the evolution Einstein equation, 
in the space of increasing and continuous functions on / = [to ,tQ + T]. Wc show 
that a has an exponential growth and the assumption of the first step on a is 
then satisfied. 

In a third step, wc consider the coupled Einstein-Boltzmann system, and, 
relying on the above results, we prove a global in time existence theorem of the 
solution (a, f) of the Cauchy problem, by applying the fixed point theorem in 
an appropriate fimction space. 

Our method preserves the physical nature of the problem that imposes to 
the distribution function f to be a non-negative function and nowhere, we had 
to require smallness assumption on the initial data, which can, consequently, be 
taken arbitrarily large. 

The paper is organized as follows: 
In section 2, we introduce the Einstein and Boltzmann equations on a Robertson- 
Walker space-time. 

In section 3, we study the Boltzmann equation in f. 

In section 4 we study the Einstein equation in a. 

In section 5, we prove the local existence theorem for the coupled Einstein- 
Boltzmann system. 

In section 6, we prove the global existence theorem for the coupled Einstein- 
Boltzmann system. 

2 The Boltzmann equation and the Einstein equa- 
tions on a Robertson- Walker space-time. 

2.1 Notations and function spaces 

A greek index varies from to 3 and a Latin index from 1 to 3, unless otherwise 

3 

specified. We adopt the Einstein summation convention aab" = ^ Uab"'. We 

consider the flat Robertson- Walker space-time denoted {M^, g) where, for 

X = (x") = {x°,x^) e M^, x^ = t denotes the time and x = (a;') the space, g 

stands for the metric tensor with signature (-, -|-, -|-, -|-) that can be written: 

g = -dt^ + a^{t)[{dx^f + {dx'^f + [dx^f] (2.1) 

in which a is a strictly positive function of t, called the cosmological expansion 
factor. We consider the coUisional evolution of a kind of uncharged massive 
relativistic particles in the time oriented space-time (]R^, g). The particles are 
statistically described by their distribution function we denote /, which is a 
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non- negative real- valued function of both the position (a;") and the 4-momentum 
P — (p") of the particles, and that coordinatize the tangent bundle r(R'*) i.e: 

/ : r(M4) ~ R4 X R-* ^ R+ , i-> f{x°',p") e R+ 

For p — (p*), q — {p^) e R'^, we set, as usual: 

3 3 

M = IpI = [E(p^)']^ (2-2) 

1=1 1=1 

We suppose the rest mass m > of the particles normalized to unity, i.e 
we take m = 1. The relativistic particles are then required to move on the future 
sheet of the mass-shell whose equation is g{p,p) = — 1. 
From this, we deduce, using H2.1|l and H2.2|l : 

/ = yrr^W (2.3) 

where the choice of p'' > symbolizes the fact that the particles eject towards the 
future . (|2.^^|l shows that in fact, f is defined on the 7-dimensional subbundle of 
T(R^) coordinalized by (a;"), (p*). Now, we consider the spatially homogeneous 
case which means that f depends only on t and p = (p^). The framework we 
will refer to will be the subspace of L-'^(R^), denote Ll(M}^) and defined by : 

Ll{R') = {f e L\M:'),\\f\\ := f vTT^|/(p)|#<+oo} (2.4) 

where \p\ is given by H2.2|l : ||.|| is a norm on L2(I^"^) ^-^^d (L2(R"^), ||.||) is a Banach 
space. 

Let r be an arbitrary strictly positive real number. We set: 

Xr = {fe Ll{R'), f > a.e, ||/|| < r} (2.5) 

Endowed with the metric induced by the norm ||.||, Xr is a complete and con- 
nected metric subspace of (L2(R''), ||.||). Let I be a real interval. Set: 

C[/; L^(R^)] = {/ : / ^ LI{R^), f continuous and bounded } 

endowed with the norm: 

Ill/Ill = ^HI/WII (2.6) 
tei 

C[I; L2{M.^^)] is a Banach space. Xr being defined by (|2.5|l . We set: 

C[I;Xr] = {/ e C[I; Ll{R% f{t) e X,. , V< G /} (2.7) 

Endowed with the metric induced by the norm |||.||| defined by H2.6() . C[/;X,.] 
is a complete metric subspace of (C[/; L2(R'^)], |||.|||). 
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2.2 The Boltzmann equation in (R^, g) 

In its general form, the Boltzmann equation on the curved space-time (R'', g) 
can be written: 

Lxf^QifJ) (2.8) 

where 

Lxf ■.^p''^-Vl,p^'p''4^ (2.9) 

denotes the Lie derivative of / in the direction of the vector field X tangent 
to the trajectories of the particles in TM^, and whose local coordinates are 
(j3", —V^^p^p^) in which F"^ are the Christoffel symbols of g; Q is a non-linear 
integral operator called the " Collision Operator" . We specify this operator in 
detail in next section. Now, since f depends only on t and (p*), (|2.8|l can be 
written using H2.9|l : 

-r»'^0 = Q(/,/) (2.10) 

We now express the Christoffel symbols defined by: 

Ta/J - ^5^''[5a5p/3 + 5/35aA. " ^P^a/s] (2.11) 

in which, the metric g is defined by (|2.1|l and g^^ denotes the inverse matrix of 
5am; iEU gives: 

.9™=.9oo = -l; =a2; 5" = a-2; = .g°' = 0; 5^ = g'^' = /or i 7^ j 



A direct computation, using H2.11(l then gives, with a 



da . 
dt ■ 



rO=aa; T^q = T^, = -; = 0; ^ for a ^ P; = 

(2.12) 

The Boltzmann equation (|2.10|l then writes, using (|2.12ll : 

in which p^ is given by (|2.3() . H2.13|l is a non- linear p.d.e in f we study in next 
section. 

2.3 The Einstein Equations and the coupled system 

We consider the Einstein equations with a cosmological constant A and that 
can be written: 

- 2^f"^ + ^^"/^ = 87rGT„^ (2.14) 

in which: 

Rap, is the Ricci tensor of g, contracted of the curvature tensor of g; 
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R = g'^f^Rajs = R" is the scalar curvature; 

Tap is the stress-matter tensor that represents the matter contents, and that is 
generated by the distribution function f of the particles by: 

j.ap ^ f tp^f^midp^dp'dp^ (2.15) 

in which l^l is the determinant of g, we have, using 12.1|l . \g\i ~ a'^. Recall that 
f is a function of t and p — (p*); then T"^ is a function of t. 
G is the universal gravitational constant. We take G = 1. The contraction of the 
Bianchi identities gives the identities S/aS"^ — 0, where Sap = Rap — \Rgap 
is the Einstein tensor. The Einstein equation H2.14|l then implies, since = 0, 
that the stress-matter tensor T"^ must satisfy the four relations VqT"^ = 
called the conservations laws. But it is proved in [S] that these laws are satisfied 
for all solutions f of the Boltzmann equation. Now if R^ ^^np ^'^s the components 
of the curvature tensor of g, we have: 

{Rap 
where 
R^ t.,aP 

in which F^^ is defined by H2.11(l . This shows that the Einstein equations H2.14|l 
are a system of non-linear second order p.d.e in g^p- In order to have things 
fresh in mind, when we will study in details the Einstein equations, we leave 
the expression of (|2.14l) in term of a, to section 4 which is devoted to this study. 
The Einstein-Boltzmann system is then the system (|2.13|l - (|2.14|l in (a, /) 

3 Existence Theorem for the Boltzmann Equa- 
tion 

In this section, we suppose that the cosmological expansion factor a is given, and 
we prove an existence theorem for the initial value problem for the Boltzmann 
equation (|2.13l) . on every bounded interval I — [to, tg+T] with to E T G 
We begin by specifying the collision operator Q in (|2.13() 

3.1 The Collision Operator 

In the instanteneous, binary and elastic scheme due to Lichnerowicz and Chernikov, 
we consider, at a given position (t, x), only 2 particles collide instanteneously 
without destroying each other, the collision affecting only the momenta of the 
2 particles that change after the collision, only the sum of the 2 momenta being 
preserved, following the scheme: 



R' 



a,\p 



^ i^Oi^ n.f3 ^ u(3^ fid 



(2.16) 



7 



{t,x) 



p + q^p' + q' 

The collision operator Q is then defined, using functions f, g on M.^ by: 

Q{f,9) = Q+{f,9)-Q-{f,9) (3.1) 

where 

Q+if,9){p)^ f ^[ f{p')9{q')A{o^,P.<i.P',q')di^ (3-2) 

Q-if,9){p)= f f f{p)9m{a,P,q,p',q')duj (3.3) 

whose elements we now introduce step by step, specifying properties and hy- 
potheses: 

1) 5^ is the unit sphere of M.^ whose volume element is denoted dw. 

2) A is a non-negative real-valued regular function of all its arguments, called 
the collision kernel or the cross-section of the collisions, on which 
we require the following boundedness, symmetry and Lipschitz continuity 
assumptions: 

< Aia, p, q, p', q') < (3.4) 

Ma,P, q, p', q') = A{a, q, p, q', p') (3.5) 

-4(a,p, q, p', q') = A{a, p' , q' , p, q) (3.6) 

\A{ai,p, q, p', q') - A{a2, p, q, p', q')\ < j\ai - aa] (3.7) 
where Ci and 7 are strictly positive constants. 

3) The conservation law p + q = p' + q' splits into: 

P° + g°=p'° + g'" (3.8) 

p + q^p' + q' (3.9) 
and (|3.8(l shows, using H2.3|l . the conservation of the quantity: 



= Vl + a' + Vl + a' I 9p (3.10) 
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called elementary energy of the unit rest mass partcles; we can interpret 
(|3.9|l by setting, following Glassey, R, T., in [?]: 



b-'=p- + 5(p-g-c.)a; ^^^^^^ 
[q' = q- b{p,q,uj)uj; iv e 

in which b{p, q, u) is a real- valued function. We prove, by a direct cal- 
culation, using H2.3|l to express q'^ in term of p' , p', and now H3.11|l 
to express p' , q' in term of p, g, that equation (|3.8|> leads to a quadratic 
equation in b, that solves to give: 

hi-- \ 2p°q°ea^uj.{q~p) 

b{p, q, uj) = — t-^-\\2 (3-12) 

— a^(uj.(p + q)y 

in which p = , q = ^ , and e is given by ((3.10(1 . Another direct compu- 
tation shows, using the classical properties of the determinants, that the 
Jacobian of the change of variables (p, q) — > (p', q') in x K.-^, defined by 
((3.11() is given by: 

'-P^ = (3.13) 

d{p,q) p°q° 

1(3.13(1 shows, using once more (12.3(1 and the implicit function theorem, that the 
change of variable ((3.11(1 is invertible and also allows to compute p, q in term 
Oip'.q'. 

Finally, formulae 1(2.3(1 and ((3.11() show that the functions to integrate in ((3.2(1 
and 1(3.3(1 completely express in terms of p,q,Lu; the integration with respect 
to q and oj leave functions Q'^{f,g) and Q^{f,g) of the single variable p. In 
practice, we will consider functions / on R x R'^, that induce, for t fixed in M, 
functions f(t) on ]R^ defined by fit){p) = f{t,p). 

Remark 3.1 1) Formulae J^S.l^) and J^H.l!^) are generalizations to the case 
of the Robertson-Walker space-time, of analogous formulae established by 
Glassey, R.T, in in the case of the Minkowski space-time, to which 
the Robertson- Walker space-time reduces, when we take a{t) = 1 in the 
metric 1^2.1}) . 

2) The expression b{p, q,uj) = uj.{p + q) used by the autor in 11,'^ and \i4Jl is 
valid only in the non-relativistic case. 

2 I - |2 I -|2 I 12 I I 2 

3) A — ke^"" ^'P' I "I'' ' , k > 0, is a simple example of functions 
satisfying assumptions ^3-4}) , i(.^.5() . 6() and ?[ )■ 

3.2 Resolution of the Boltzmann equation 

We consider the Boltzmann equation on [to, to + T] with to E M+, T E and 
a is supposed to be given and defined on [to, to -\-T]. 
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The Boltzmann equation H2.13|l is a first order p.d.e and its resolution is equiv- 
alent to the resolution of the associated characteristic system, which can be 
written, taking t as parameter: 

We solve the initial value problem on / = [<o, + T] with initial data: 

p\t^) = z/; /(io) = /to (3.15) 
The equation in p = (p*) solve directly to give, setting y = (y*) e M"^; 

P{h + t,y)^^p^y, ie[0,T] (3.16) 

The initial value problem for f is equivalent to the following integral equation 
in f, in which p stands this time for any independent variable in IS?: 

pto+t ^ 

f{to+t,p) = ftM+ -Q{fJ)is,P)ds te[0,T] (3.17) 

Jto P 

Finally, solving the Boltzmann equation H2.13|l is equivalent to solving the inte- 
gral equation (|3.17|) . We prove: 

Theorem 3.1 Let a be a strictly positive continuous function such thata{t) > | 
whenever a is defined . 

Let ft„ G L2(^^); /to — ^-^j ^ € Rq such that r > ||/to||. Then, the initial 
value problem for the Boltzmann equation on [to, to + T], with initial data ft^, 
has a unique solution f £ C[[io, ^o+2^]; ^r]- Moreover, f satisfies the estimation: 

Sup ||/(t)|| < ll/tJI (3.18) 

t<^[to,to+T] 

Theorem 3.1 is a direct consequence of the following result: 
Proposition 3.1 Assume hypotheses of theorem 3.1 on: a, ftg and r. 

1) There exists an integer no{r) such that, for every integer n > no{r) and 
for every v £ X^, the equation 

\/nu TrQ(u,u) = v (3.19) 

po 

has a unique solution Un G Xr 

2) Let n eN, n> no{r) 

i) For every u € X^ , define R(n, Q)u to be the unique element of X^ 
such that: 

y/nR(n, Q)u - \q [R{n, Q)u, R{n, Q)u] = u (3.20) 
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(3.21) 



a) Define operator Qn on Xr by: 

Qn{u, u) = n^JnRin, Q)u — nu 
Then 

a) The integral equation 

fito + t,p) - ftM + / Qn{fJ){s,p)ds t e [o,r] 

Jto 



has a unique solution fn G C[[io,^o 
Moreover, /„ satisfies the estimation: 

Sup < ll/toll 

t£[to,to+T] 



(3.22) 



(3.23) 



b) The sequence {fn) converges in C[[to,to + T]; Xr]to an element 
f G C[[to,tQ+T]; Xr], which is the unique solution of the integral 
equation p. j7| ). The solution f satisfies the estimation ^y.iyp . 

Proof of the proposition 3.1 

The proof follows the same lines as the proof of theorem 4.1 in JS]. We will 
emphasize only on points where differences arise with the present case and we 
show how we proceed in such cases. 
Proof of point 1) of prop 3.1 
We use: 

Lemma 3.1 Let f, g e Ll{R^). then j,Q+ {f , g), ^Q- {f , g) L\{R^) and 



<C{t) II /nil 5 



-oQ"(/,5) 



<C(i)||/||||g|| (3.24) 



i^Q^ifJ) - 4TQ+(ff,5) < cm / II + II 3 II) II .f - 9 II 



^Q-{fJ) - ^Q-{9,9) < C{t)i\\ / II + II .g II) II /-5 II 



^QifJ)~^Q{9,9) 



(3.25) 

(3.26) 
(3.27) 

Vl+W < \/l + < \/l + + = \/l + ^{p" + q'") < 2(p'" + g'°) 



where 



<C(t)(||/|| + || 5 II) 11/ -5 II 
C{t) =32TTCia^{t) 



proof of lemma 3.1 

We deduce from H3.8|l and a > 1 that: 



1 
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Expression (|3.2|l of Q'^{f,g) then gives, using (|3.4|l : 

dp< 2a^{t)Ci 



if, 9) 



— -0 — Q if^a) 



\f{p')\\g{q')\dpdqdcj 



We then deduce, using the change of variables H3.11|l and H3.13|l that gives 
dpdq 



^dp'dq' and the fact that, by ^ 



< 2 



The estimation of 



\f{p')\\g{q')\dp'dq'<C{t)\\f\\g\\ 



^Q-{f,9) 



follows the same way without change of vari- 
ables and H3.24|l follows. The inequalities H3.25|l are consequences of (|3.24|) and 
the bilinearity of and Q~, that allows us to write, P standing for or 



P{f. I) - P{9: 9) = P{f: f~9) + P{f~9,9)- 

Finally, (|3.2()|l is a consequence of (|3.25|) and Q = — Q . This completes 
the proof of the lemma 3.1. 

Now the continuous and strictly positive function t a?'{t) is bounded from 
above, on the hne segment [to,^o + T], and C(t) given by (|3.27|l is bounded 
from above by a constant C(to, T) > 0. Hence, if we replace C(t) by C(<o, T) in 
the inequalities in lemma 3.1, we obtain the same inequalities with an absolute 
constant than the inequalities in proposition 3.1 in |15j . The proof of the point 
1) of prop 3.1 is then exactly the same as the proof of prop 3.2 in |16j. ■ 
Proof of the point 2) of prop 3.1 
We use, no(^) being the integer introduced in point 1: 

Lemma 3.2 We have, for every integer n > no(r) and for 

\\ ^/PiR{n,Q)u\\ = \\u\\ (3.28) 



Proof of lemma 3.2 

H3.28|l is a consequence of: 



Q{fJ)ip)dp = 0, V/eii(R3) 



(3.29) 



we now establish. It is here that assumption (|3.5|1 . (|3.()|l on the collision kernel 
A are required. Define operator Q* , for f,g& Ll(R^) by: 



'{f,9){p)^l 



a^A{a,p,q,p',q') 



R3 Js^ Q 
fiP)9iq) - f{q)9{p)]dqduj 



[fip)9iq) + f{q)9{p) 



(a) 



Clearly, Q*{f,f) = Q{f,f), where Q is the collision operator defined by 1)3. 1|) . 
H3.2|l . (|3.3() . Now let $ be a regular function on M^, such that: 



^p) + a>(g) = $(p') + q,{q') 



(b) 
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Multiplying (a) by and integrating on K'^ yields: 

- f{p)9{q) - f{q)g{p)]<^{p)dpdqdu 

(c) 

Let us compute the integral of the r.h.s of (c) using the change of variables 
(j),q) ^ {<j,p)', formulae (|3.12|) shows, since p^q^ remains unchanged, that 
b{p, q, uj) change to —b{p, q, u) and formulae (|3.11() then show that {p' , q') change 
to {q',p')', next assumption H3.5|l on A show that A{a,p,q,p' ,q') remains un- 
changed. We then have, since the Jacobian of the change of variables {p, q) 
{q,p) is 1 

^Q*if,g)<i>iP)dp^^ [ M^'Pff' + fiq')gip') 

- f{p)g{Q) - f{q)9{P)]Hq)dpdqduj 

(d) 

The sum of (c) and (d) gives: 

lQ*(/,g)$(p)dp- = ^ / 5(«,p-g-f'g") [/(p-0g(g-0 + f{q')g{p') 

R3 4 Jr3xR3xS2 P T 

- f{p)9{q) - f{q)g{p)]dpdqduj 

where 



(e) 



S{a,p,q,p',q')^A{a,p,q,p',q')Mp) + ^q)] (f) 
we can write using, (e) 

/ ^Q*{f,g){p)mdp = Io-Jo (g) 

JK3 P 

where 

'^0 = ^ /r3 /r3 /s. + /(g-O.gb-Old^'- 

Jo^i /r3 /r3 /s. + fiq)g{P)]dpdqdu; 

Let us make in Iq the change of variable {p,q) {p',q') defined by H3.11|l 
and whose Jacobian (|3.13|) gives dpdq = ^o'^/o dp'dq' ; we have, since (p, is 
expressible in term of {p' , q') by virtue of H3.13|l : 



h = - I d^ I ^ ' ,0 [/(P0g(g0 + /(9').9b-')]#'rf9' (h) 



^ f f S{a,p,q,p ',q') ^ 

2 Js^ yR3xR3 J3'°9" 

Now, assumption H3.6() on A and assumption (b) on $ imply that S defined 
by (f) satisfies the relation: S{a,p,q,p' ,q') — S{a,p' ,q' ,p,q). Replacing in /q 
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S{a,p, q,p' , q') by S{a,p' , q\p, q) , we deduce, from the expression (h) of /q that 
Iq = Jo- then (g) imphcs: 



/ ^Q*{f,gMp)dp = 



(i) 



Now the function — p^ — -^/l + a^jpp satisfies hypothesis (b) as a 

consequence of the conservation law H3.8|l . The relation H3.29(l then follows from 
the above choice of ^{p), (i) and the relation Q*{f, /) = Q{f, /)■ 
Now let us prove 



We have, multiplying equation (|3.20() by = ^/l + a^\p\^, integrating over M"^, 
and using H3.29|l : 



^/l + a^\p\'^Rin,Q)u{p)dp ^ / y'T+a^uip)dp (3.30) 



If we make in each side of (|3.30ll the change of variables q — Bp where 

B = Diag{a,a,a), then |gp = ; dp — ^dq and (|3.3U|) gives, using 

definition EH of 1 1. II : 

^/^\\R{n,Q)uoB-^\\ = \\uoB-^\\ (3.31) 
But, if we compute ||uo_B||, using the above change of variable, we have: 



uoB\\= i yjl + \p\'^uoB{p)dp = 



l + -\q\2\u{q)\-dq 



1 

< — 



1 + ^ / ^l + \q\^u{q)dq 



i.e 



\uoB\\ < i- ( 1 + i 



The assumption a > | implies ^ (l + ^) < 1, so that ||uoi?|| < ||u||, this 
implies that uoB £ Xr if m G Xr- Now since (|3.31|) holds for every u G Xr, we 
have, replacing in (13.31(1 u by uoB, \\y/riR{n,Q)u\\ = \\u\\. We then have H3.28|l 
and lemma 3.2 is proved. 

Now (I3.28|l is exactly equality H3.10|l in proposition 3.3 in ^H], we then prove 
exactly as for prop. 3.3 in that all the other relations of that proposition 
hold in the present case. Using this result, the proof of point 2)a) of proposition 
3.1 is the same as the proof of prop 4.1 in [16j and the proof of point 2) b) of 
prop 3.1 is the same as the proof of theorem 4.1 in ^H), just replacing, [0, +oo[ 
by [to , to + T] . This completes the proof of prop 3 . 1 which give directly theorem 
3.1. ■ 
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4 Existence Theorem for the Einstein Equations 

4.1 Expression and Reduction of the Einstein Equations 

We express the Einstein equations H2.14|l in terms of the cosmological expansion 
factor a which is the only unknown. We have to compute the Ricci tensor Rap 
given by f lTT^ . 

The expression H2.12|l of F^^ shows that the only non-zero components of the 
Ricci tensor are the Raa and that Rn = R22 = Rss- Then, it will be enough to 
compute Rqo = R^ q ;^Qand Rn = R'^ i^^. The expression ( 12.12(1 of F^^ and 
formulae (|2.1ti|) give: 



0; ^0, 10 — -^0, 20 — -^0, 30 



^1,01-""; ^1,11-0; i??,2i - ^1,31 - (a)^ 
We then deduce that: 



00 



-3— and Ru — ad + 2{d) 



We can then compute the scalar curvature R to be 

R = R2^ g"^Ro.p = 6 



The Einstein equations ((2.14|) then take the reduced form, using expression ((2.1|l 
of g: 



too 



-{df -2ad + a'^A = SttTh 
that can be written, using r"'^ = .g^^y'^TA^: 



a 3 3 



(4.1) 
(4.2) 



d 47r .„nn 9„iix A 

- = (T"° + 3a2r") + - 

a 3 3 

in which T"^ is defined in term of f by 12.15|l . In this paragraph, we suppose f 
fixed in C[[0,T];X,], with /(O) = /o G Ll{M.^), /o > a.e. and r > ||/o||. 



4.2 Compatibility 

The relations Rqi ~ 0, Rij = ii i ^ j, Rn = R22 = R33, imply for the Einstein 
tensor Sap — Rap — ^9apR, that: 



Til — Tx 



Toi = 0, Ty =0 for j 



(4.3) 
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But the stress- matter tensor Tap is defined by (|2.15|l in terms of the distribution 
function f. So, the relation (|4.3(l are in fact conditions to impose to f. we prove: 

Proposition 4.1 Let ft„ and r > be defined as in theorem 3.1. Assume that, 
in addition, ft„ is invariant by S03 and that, the collision kernel A satisfies 

A{a{t), Mp, Mq, Mp', Mq') = g,p', q') , VA/ e ^03 (4-4) 

then 

1 ) The solution f of the integral equation satisfies: 

fito + t,Mp)^ f{to + t,p), Vte[0,T], VpeR^ VAf e (4.5) 

2) The stress-matter tensor Tap satisfies the conditions \4.3\j . 
Proof 

1) Let M e S'os; *^n7ll gives, since ft^Mp) = /(^(p) and / = 

f{t^ + t,Mp) = ftM + -o^QifJKs,Mp)ds,te [0,T] (4.6) 

Jto P oM 

Notice that p° — p° (p) is invariant by Soa ■ Now, definition , H3.2|l , H3.3|l of 
Q gives: 

Q{f,f){s){Mp)= f ^li^ f [f{s,p')f{s,q')-f{s,Mp)f{s,q)]A{a,Mp,q,p',q')duj 

(4.7) 

Let us set in (|4.7|) q = Mqi; lo = Muoi. Then formulae H3.11|l give using 
expression (|3.12|l of b, the invariance of the scalar product in by 6*03 : 

{ p' = Mp + b{Mp,q,uj)ui = Mp + b{Mp, Mqi, Muji)Muii = Mp + b{p,qi,uji)MuJi 
\ q' = q- b{Mp, q, uj)uj ^ Mqi~ b{Mp, Mqi, MuJi)Mu;i ^ Mqi + b{p, qi,uJi)Muji 

So that p' = Mp'i, q' = Mq'i where: 

Pi =p + b{p,qi,uji)uji; q[ ^ qi - b{p,qi,uJi)uJi 

Then , (|4.7|) implies, using assumption (|4.4|) on A, = gj, and the invariance 
of the volume elements dq, dto by 6*03 

Q(/, f)(s){Mp) = Q[f{s)oM, f{s)oM]{p) 

H4.6|l then write: 

/(io + t)oM{p) = ftM + ^Q{f(s)oM, f{s)oM){p)ds 

Jto P 

which shows, by setting h{s) — f{s)oM, that ||/i(s)|| = ||/(s)|| setting q ~ Alp 
in the integral in p defining ||/(s)oM|| , and that, h and f are 2 solutions in 
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C[[to,to + T];Xr] of the integral equation (|3.17l) : the uniqueness theorem 3.1, 

then imphes h = f and the first point in prop 4.1 is proved. 

2) Let us consider the foUowing element of S03 in which 6* e M^: 



Mi(0) = 




Consider the expression (|2.15|) of T"^ in which f satisfies (|4.5I) and observe that 
P°[Mk{9)q] = k = 1, 2, 3, e M. 

(i) Set in (|2.15() a = 0, P — i with i = 1, 2; now compute the integral using 
the change of variable p — Mi(7r)g; the integral in q gives, using H4.5|l : 

= ~T°\ i = 1, 2; hence = = 0. 

(ii) Set in (|2.15l) a = 0, /3 — 3 and compute the integral in q using the change 
of variable p — M^{7T)q; the integral in q gives, using 1)4. 5|l : T^^ = —T^^ 
hence ^0. 

(iii) Set in H2.15|l a — i,l3 = ji^j and compute the integral using the change 
of variable p = Mk{f)q with k = 1 if i = 1, j = 2; k = 2 if i = 1, j = 
3, k = 3 if i = 2, j = 3 the integrals in q give, using T^^ ^ _yi2. 

yl3 ^ _yl3. ^ _j.23 i^g^^g ^ yl3 ^ y23 ^ q. 

(iv) Set in (|2.15l) a = (3 — i with i = 1, 2; and compute the integral using the 
change of variable p = Mk{^)q; taking k = 1 if i = 1 and k = 3 if i = 2, 
to obtain using = T^^ and T^^ = T^^. 

This completes the proof of prop 4.1. 

In all what follows, we assume that ftg is invariant by Soa and that the 
collision kernel A satisfies assumption 1)4. 4|l . Notice that A defined in Remark 
3.1 is an example of such a kernel. 



4.3 The Constraint Equation 

We study the Cauchy problem for the system H4.1(l - I|4.2|l on [0,T] with initial 
data: 

a(0) = ao; a(0) = do (4.8) 
The Einstein equations (|2.14() with G = 1, give, using the Einstein tensor S"^: 

Hi := Si + A.gO - 8^T^ = 

it is proved in P- 39, that, in the general case, the four quantities Hi do 
not involve the second derivative of the metric tensor g with respect to t, and, 
using the identities Va(S'"^ - Ay"'' - SttT"'') = 0, that the four quantities Hi 
satisfy a linear homogeneous first order differential system, consequently: 

1) For < = 0, the quantities Hi express uniquely in term of the initial data 
ao, do and /o. 
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2) If iy°(0) = 0, then H^(t) = in the whole existence domain of the solution 
g of the Einstein equations. 

But -ff^(O) — are 4 constraints to impose to the initial data at i = and 
the equation 7?° = are called constraint equations. In our case, the rela- 
tions 5° ~ 0, r" = (see prop 4.1) show that the constraints — 0, called 
the momentum constraints are automatically satisfied. It then remains the 
constraint with a — which is equivalent to S*™ — A — 87rT°° = 0, called the 
Hamiltonian constraint, and which is exactly H4.1|l . 

So, following what we said above, the Hamiltonian constraint (|4.1() is satisfied 
in the whole existence domain of the solution a of the initial values problem on 
[0,r], if the initial data aq, dp, /o at t = 0, satisfy, using expression (|2.15|l of 
r"*', the constraint 

H4.9|l gives two possible choices of do, when ag and /o are given. We will choose, 
taking also into account the hypothesis on a(t) in theorem 3.1: 

ao > I /o e i^(R^); /o > a.e do > 0. (4.10) 
We now concentrate on H4.2I) which is the evolution equation. 

4.4 The Evolution Equation 

We set e = 3f , then ^ = 3[f - (^f] and lO gives: 

^ = -— -47r(T°° + 3a2r")+A (4.11) 
3 

H4.11|l is the Raychaudhuri equation in 6. We prove: 

Proposition 4.2 There can exist no global regular solution for the coupled 
Einstein- Boltzmann system in the case A < 

For the proof, we use the following result, proved in 1^: 

Lemma 4.1 Let u and 6 be 2 differentiate functions of t satisfying: 

^<-y; " = -y; e{t^)^u{t^) (4.12) 

for a given value ti of t. Then 0{t) < u{t) for t > ti. 
Proof of proposition 4.2 

Suppose A < 0. Then the Raychaudhuri equation (|4.11(l gives 9 < We 
have, integrating the equation in u on [ii, t] when u ^ and since 0(ti) — u{ti): 
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Let us show that we have necessary a < 0. 



The Hamihonian constraint 14.1(1 writes, using and p° = v/l + a^(i)|pp- 



The derivative of the l.h.s is: A = 2(^) ii - (S-) 



the derivative of the r.h.s is 



Sa^a/ p"f(t,p)dp + a^ [ ^f{t,p)dp+f p'%{t,p)dp 

Jr3 P^ Jr3 dt 



Bit) = ^ 
^ ' 3 

(b) 

But we have by ipmjl ^ = ^^j^, then IpO?^ shows that the last integral in 
B(t) is zero. We have, naturally, A(t) = B(t). Suppose d{ti) > 0; (b) imphes 
B{ti) > 0. Now ^ < ^ ^ < 0. But 6* = 3f. Then ^ = 3 [f - < 0, 

so that A{ti) < and the equality A{ti) = B{ti) is impossible. We then have 
d{ti) < 0. The hypothesis d{ti) — would implies that the r.h.s of (a) be a 
constant, but it is not the case, since it depends on f which changes with ftg 
We then conclude that, necessarily a < 0, which implies < 0, so that in H4.12|l 
we have 9{ti) < 0. Now (I4.13|l shows that , since by ((4.12|l u is a decreasing 
function: 

3 

u{t) —oo when 3 + 9{ti){t — ti) 0, and ti < t < ti — — — r- := t* 

> 0[ti) 

(4.14) 

By lemma 4.1, (|4.14ll implies that 9{t) = 3^(i) — > -oo when t t*, then: 

+00 when t t* (4.15) 
ait) J < 

Now, since a < 0, a is a decreasing function on and then, a{t) < a{ti), 

yt g (a) and l|4.1|l then implie, since A < 0: 

(af < Y^Hti)V^ + aHh)\\\f\\\ := C\hJ) , Vt E [ti,t*[ 



so that: 



W)) 



H4.15|l then implies that ^2 (V) ^ when t ^ t* and this can happen only 

if: a{t) ^ when t ^ t* . So the cosmological expansion factor a tends to zero 

in a finite time and such a solution (a, f) cannot be global towards the future. 
This complete the proof of proposition 4.2. ■ 

We now study the case A > 0. Notice that since r"° > 0, > (HlJ-g^ 
gives by subtraction ^ — {^)'^ < 0; this implies ^(f ) < 0, which shows that, in 
all the cases, = 3- is a decreasing function. In the case A > 0, the Hamiltonian 
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constraint (|1?T|) gives: 
equivalent to: 



>f i.e 



a(t) 
a{t) 



ait) I / A^ 
a{t) V 3) 



a{t) - V 3 

a(t) /A 
a(i) - V 3 



> 0, which is 



(4.16) 



(4.17) 



The continuity of i — > imphes that we have to choose between (|4.16f) and 
H4.17|l . Since a > 0, 1)4. 17|) imphes a < and a is decreasing; H4.16|l imphes that 
d > 0, then a is increasing and since ^ is decreasing, this gives on [to,t]: 



a{t) > a{to); 



3 ~ a{t) ~ a(to)' 



t > 







(4.18) 



RecaU that our aim is to study the coupled Einstein-Boltzmann system and 
we had to require, for the study of the Boltzmann equation, that a, which is 
positive be bounded away from zero. This problem is solved by choosing (|4.16() . 
Another important consequence of (|4.16|l is that it implies on [to,t] 



a(t) > a(io)eV 3(*-*o), t > to 

Which shows that, the cosmological expansion factor has an exponential growth; 
but, it also shows that, an eventual global solution a will be unbounded, and 
this is why, in order to use standard results, we make the change of variable: 



which gives: 



(4.19) 



(4.20) 



Then, using lOsl T™ and T^^ in the r.h.s of 14.21 express in, terms of e and f, 

and we set: 

1 



00 



1 



1 



\p\\f{t,p)dp 



P = a^T 



2rpll 



1 



dp 



+ 

p stands for the density and P for the pressure. One verifies that: 

P <p 

Recall that r > is such that r > ||/o||. If we set: 



do 



167r 



(4.21) 
(4.22) 

(4.23) 
(4.24) 
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One checks easily, using a > ao > |, (|4.18|) with to = H4.10|l H4.9|l and (|4.24|) . 
that 

e^i-GlO,-]; et 6* = 3- G [VSA , do] (4.25) 

a 3 a 

A direct calculation shows, using ^J^, (jOUjl . l|OT|l . llO^ , and KTT]i . that 
the Einstein evolution equation H4.2|l is equivalent to the following first order 
system in {e,6): 

Q 

e = --xe (4.26) 
o 

n2 

e = -—-47:{p + 3P)+A (4.27) 
with p = p{e, /) and P P(e, /) given by KT^i and 

we will study the initial values problem for the system (|4.2()f) - (l4.27l) with initial 
data (60,6*0) and t = 0. By virtue of the change of variables (|4.25() . we will 
deduce solution for the Einstein evolution equation H4.2|l by setting: 

e(0) = — e{0) = 3^ (4.28) 
ao ao 

in which oo, ao satisfy the constraint H4.10|l - H4.9|l with /o given. By virtue of 
H4.25|l it will be enough for the study of the initial value problem for H4.26|l - 
H4.27|l to take cq, ^o such that 

< eo < I < 00 < rfo (4.29) 



Remark 4.1 The equivalence of the evolution equation i4.S\ ) and the system 
14.26]) - ^~FJ\ I requires that any solution (e, 0) of l4.26]) - ^~F^\ )- i4.29i l satisfies 
e > and 9 > 0. 

In fact, l|OT|) . require e ^ 0; by g^HI), e(0) = eo > 0. Since e is 

continuous and cannot vanish, by virtue of the mean values problem, e remains 
strictly positive. Now if we have 0{ti) ~ 0, then (|4.26() imphes e(ti) — and 
l|4.2()|l implies d(ti) = 0, which is impossible given (|4.16(l . So 9 never vanishes 
and by l^TM 9{0) = ^o > 0; So 6* > 0. 

For the global existence theorem, we will need the following a priori estimation 

Proposition 4.3 Let S > and to E be given; suppose that in 14.21}) and 
{4-22^ , f e C[to,to + S; L2(M.^)] is given. Let (e, 0) be any solution of the system 
\4.2b\ - ^4-'i'^ on [to, to + S]. Then 9 and a = \ satisfy the inequalities: 

dito + t) < 9{to) + A, te[0,S] (4.30) 
a{to + t)< a(to)e(^+*)(*«+*+i)', t e [0, S] (4.31) 
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Proof. 

Since p > 0, P > 0, H4.27() implies: ^ < A. Integrating this inequality on 
[to, to + t] where t G [0, (5] yields (|On|) . 

Next, since e > 0, integrating (|4.2t)|) that writes — f = f over [toi ^o+^]) i G [0, S\ 
gives: 

a(to + t) < a(to)e^'o ^* ^''^^ t G [0, 6] (4.32) 
Now setting in H4.3()|l s = to + i G [^o, + <5] ,and integrating on [to, to + 1] yields: 

f^'^''^ds<['^ + ^){t,+t + l)^ te[0,<5] (4.33) 
H4.31|l then follows from (a) and (b).H We deduce: 

Proposition 4.4 Let T > and f E [[0,T];Xr] be given. Suppose that the 
initial value problem \4.26}) - l4-'i^ - iJ-2S{ ) with initial data | |^.^<y| ) at t — sat- 
isfying the constraints i4.10j) - ^4-S^ ^0,^ o, solution (S = jj'®) [0,to] with 
< to < T. Then, any solution (e = -^,9) of the initial value problem for the 
system \4-2b\ - ^.21\j on [to,to+(5], 5 > 0, with initial data {e,9){to) = (S, 0)(to) 
at t ^ to, satisfy the inequalities: 

a{to +t)< C2e^^(*°+*+i)', t e [0, S] (4.34) 



^(to +t) < 371 + A(r + t), te[0,S] (4.35) 

where: 



C2 = aoe^- ^3=71+3^; 71 = 7i(ao, r, T) = (- + ^ - + 3ra^)(r + l)^ 

(4.36) 

Proof 

We apply proposition 4.3 to the solution {e,9) of H4.26|l - H4.27|l on [to, to + S]; 
l|4.31|l gives, since a(to) = S(to), 6'(to) = 6(to): 

a(to + t)<S(to)e(^+*)(*"+*+l)^ te[0,S] (a) 

Now apply l|OT|l to the solution (2,6) of l|l?^ - (|07|) on [0,to], at the point 
to- We obtain by setting in H4.31|l . to = 0, t = to and since a(to) — S(to), 
S(0) = ao, e(0) = 3f^, < to < T: 

S(to) < aoe(^+^)(*°+^)^ < aoe(^+*)(^+^)^ t e [0,S] (b) 

Now since the initial data at t = satisfy (|4.1()|l . 14.9|l . this imply, using ao > 0, 
ao>l, ||/o||<r: 




(c) 
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and (b) gives: 

S(to) < aoe^' (d) 

with 

Now, we apply (|On|) to the solution (S, 9) of (|0^ - (|07jl - l|i:^ . on [0, to] at 
point to, t = to and since S(0) = 3f^, 0<to <T: 

E{to) < S(0) + Ato < 3^ + AT (f) 
We obtain, using (c) and (f): 



^+3^(^3+V3+T™°j + 3^ 
This gives, using the definition (e) of 71: 

E(to) A A 

^+3<7i+3T (g) 
and (|4.34|) follows from (a), (d) and (g). 

The inequahty follows from l|On| . 6'(io) = S(io), (f), (c) and (e). This 

completes the proof of prop. 4.4. ■ 
We can now prove. 

Proposition 4.5 Let T > and f £ C[[0,T]; Xr] be given. Then the initial 
value problem for the system | |^.i?6p - f4.i?7| ) with initial data (eo,0o) at t — 
satisfying \4..2y\j , has an unique solution (e, 9) on [0, T] 

The proof of proposition 4.5 will use the following result. 

Lemma 4.2 Let ei, 62 g]0, |]; then we have, C being a constant: 

|p(ei,/) - p(e2, /)| < 4^||/(t)|||ei - 62! (4.37) 
6162 

|P(ei,/) -P(e2,/)| < 4^||/(t)|||ei - 62! (4.38) 
Proof of the lemma 4.2 

1) To prove, H4.37|l . we write H4.21(l in ei, 62 e]0, |] and we subtract. We will 
use: 



1 1 



'1 - 62! |ei - 62! I 2 , , 2i ^ o lei - 62 ^^ 

^3^3 = ^3^3 |e2 + eie2+ei| <3 3 3 (a) 
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We will also use, applying \\/a — \/b\ = |a — 6||-ya + Vb\ ^, < < 1 and (a): 



< 



< 



1 1 



1 + 



(b) 



3y/T+W\ei-e2\ 



Now (IOT|) gives: 

|p(ei,/)-p(e2,/)| < 



1 1 



ll + -\p\2f{t,p)dp 



'1 + -\P\^-. l + -\p\2 



(c) 



f{t,p)dp 



and (|4.37|) follows from (a), (b) (c) and < 1. 



2) To prove, 
will use: 

1 1 



we write (|4.22|) in ei, 62 €]0, |] and we subtract. We 



5^5 



eVe: 



5 5 .et+e^ea+e^e^+eie^+e^l < S^^i^^ 



1^2 



eYe: 



(d) 



and we also use < < 1 and | -/^ < 1 to obtain: 



(pi 



1^2 



< 



< 



1 



1 



ip^ 
\p\' 



1 1 



< ^ Jl + ^|p|2|ei - esl < -^^l + \p\^\ei - 62! 



Now (IOT|) gives: 

|^^(ei,/)-P(e2,/)| < 



1 1 



1 + 

(pl)2 (p2)2 
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f{t,p)dp 



(e) 



(f) 



f{t,p)dp 



H4.38|l then follows from (d),(e), (f), < < 1, and using in the first integral 

(pl)2 I ] 



i + j^\P? 



<e?Jl + 3|pP<eiyrW 



(4.39) 
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This completes the proof of the lemma 4.2. ■ 
Proof of proposition 4.5 

Define using 

gSni, (F^ . the function F on ]0, |] x]0, do] by: 



F{e,e) = 



-— -47r(p + 3F) + A 



(4.40) 



To prove the proposition 4.5, we first prove the existence of a local solution. To 
have that result, we have to show that F defined by H4.4()|l is continuous with 
respect to t, and locally Lipschitzian in {e,9) with respect to the norm of R^. 
The dependence of F on t is through f that appears in p and P (see formulae 
(ji:^ and KT^ ). But, since / S C[to,to + T,Xr], the definition (ESI of this 
space shows that F is a continuous function of t. Now we have, with 
(ei, 6*1) , (e2, 02) G ]0, l]x]0, do], and by usual factorization: 



71 ei 
3 



6*262 



< 



l + do 



ei - 62 



< 



Mo 
3 ' 



Let us fix e]0, 3] and we take 61,62 G 
and we deduce from H4.37|l - H4.38|l that 

|p(ei,/)-p(e2,/)|+3|P(ei,/) 



then 



/'(62,/)| <^(6l)r|6i 



< 



62 



(4.41) 



(4.42) 



1,2 



(4.43) 



Where N is a constant depending only on . 
Then, if is endowed with the norm ||(x,y)||R2 
defined by (|On|) . using (|0T|) - (|0^ - (|0^ : 



we have for F 



\\F{ei,0i) - F(e2, 02)||r2 < M(|^?i - + \ei - 62!) = M\\{ei,0i) ~ (61, 0i)||r2 

(4.44) 

Where M is a constant depending only on e^, and do. 

(|4.44|l shows that F is locally Lipschitzian in (e, 0) with respect to the norm 
of M^. The standard theorem of the first order differential system implies that 
the initial value problem H4.26|l - I|4.27|l - I|4.29|l has a unique local solution {e,0) 
on [0,(5], 6 > 0. Notice that since / e C[[0,r];X^], the system (jl^Hl-lIl!^ is 
defined on [0,T]. Hence the maximum value oi 6 is 6 = T. Now if < (5 < T, 
prop. 4.5 in which we set to = Q, Q<t<S<T, shows, applying H4.31|l and 
a — ^, that for every solution e, the function i is uniformly bounded, and, p 
and P defined by (|4.21|) and (|4.22|) are uniformly bounded. Consequently, the 
function F defined by H4.40|l is uniformly bounded. The standard theorem on 
the first order system then implies that the solution {e,9) is global on [0, T]. 
This complete the proof of prop. 4.5. ■ 

We deduce a result that will be useful to prove the global existence for the 
coupled Einstein-Boltzmann system. We will use the number Do defined by: 



Do = 371 + A(T 
where 71 is defined by H4.36|l . We prove: 



1) 



(4.45) 
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Proposition 4.6 Let T > and f € C[[0,T];Xr] be given. Let (S = ^,6) be 
the solution of the initial value problem for the system j^.i^6] )- f^.g7| ) with initial 
data {eo,9o) at t — given by i4.28)) with ao, clq satisfying the constraints 
i4.10\) - U^ . Let tci € [0,T]. Then the initial values problem for the system 
i4.26)) - ^4.27}) , with the initial data (S = , 0)(to) att = to has a unique solution 
(e — 6) on [to, to+5], where 6 > is independent oft^. The solution (e = ^,6) 
satisfies the inequalities: 

I < aito + t)< C2e^^(*«+*+i)' , te[0, S] (4.46) 

V3A<9{to+t) <3ji+A{T + t), te[0,S] (4.47) 
where C2, C3 and 71 are defined by 
Proof 

We have 6*0 — and H4.10|l implies 6*0 > \/3A, so, the proof given for prop. 4.5 
with the function F given by H4.40|l and defined this time on ]0, |] x [-v/SA, do] 
leads to a solution (S, 9) satisfying 9 > \/3A, hence, 9(io) > \/3A. Now 
suppose that we look for solutions (e, 9) on [to, to + 5] with < (5 < 1; then 14.35|l 
shows that every solution {e,9) satisfies 9{to + t) < Do, < t < 6, where Do is 
defined by H4.45|l . We prove the existence of a local solution (e, 9) of H4. 26(1 - 14. 27|l 
on [^0,^0 + S], with initial data at t — to', eito) = S(to), 9{to) = 9(fo) > \/3A, 
following the same lines as in the proof of prop. 4.5, with the function F given 
by (|4.4U|) . defined this time on ]0, |] x [-s/SA, Do]- This leads to the existence of 
a local solution (e, 9) on some interval [to,to + 5], Q < 5 < 1. On the other hand, 
lH:!^ gives, using < to < T , Q < 5 < I: a{to + t) = ^^^^ < C2e^^^'^+'^^\ 
which shows that, for every solution (e, 9), -i is uniformly bounded. From there, 
the existence of a number < 5 < 1, independent of to- Finally, 14.46|l follows 
from (|4.34() and ^^.^^^^^^ = e{to+t) < |; 14.4711 follows from H4.35|l and what we 
said above. This ends the proof of prop. 4.6. I 

Theorem 4.1 Let T > and f e C[[0,T];Xr] be given. Then the initial 
value problem for the Einstein equation j^. j| )- f7^ - i^-'j^ with initial data ao, do 
satisfying the constraints j^. has an unique solution a on [0,T]. 

Proof 

Apply proposition 4.5, choosing the initial data (eo, 6*0) at t = by (|4.28() : H4.25|l 
shows that (60,6*0) satisfy H4.29|l . Prop. 4.5 then proves the existence of the 
unique solution {e,9) of (|0^ - (|07I) - (|0S1) on [0, T]. Now the equivalence of 
the evolution equation H4.2II with the system (|4.26|) - (|4.27|l show that a = -i is 
the unique solution of (|4.2|) with the initial data ao, ao- Now since ao, ao sat- 
isfy the initial constraint (|4.10|l . The Hamiltonian constraint H4.1() is satisfied 
on the whole existence domain [0, T] of a. This ends the proof of theorem 4.1. ■. 
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5 Local Existence for the Coupled Einstein-Boltzmann 
System 

5.1 Equations and Functional Framework 

Given the equivalence of the initial value problems H4.1|I - 14.2|) - (|4.8I) with con- 
straint (|ni7|) - (|OI) and (|0^ - (|07I) - (|0S1) . and following the study of the Bohz- 
mann equation , paragraph 3, by the characteristic method that leads to (|3.14|) . 
we study the initial value problem for the first order system: 



dt 



-i: = -oQifJ) (5-1) 



e = --e (5.2) 

^ = -y-4^(p + 3P) (5.3) 

with, in P = p{ej) and P = P{eJ) defined by (|OT|) and The 

initial data at i = are denoted /o ; eo ; i.e. 

/(O,p) = /o(p); e(0) = eo; 0(0) = (5.4) 

where we take 

foGLliR^); /o>0 a.e;eo = — ; 0o - 3^ (5.5) 

ao ao 

With, following H4.10|l - H4.9|l . ag, ao, fo subject to the constraints: 



ao>l - = J T / «^\/l + «^bP/o(p)# + ^ (5.6) 

We are going to prove the existence of the solution (/, e, 6) of the above initial 
value problem, on an interval [0, 1], I > 0. Given the study in paragraphs 3 and 
4, the functional framework will be the Banach space E = L2{M.'^) x R x R, 
endowed with the norm 

||(/,e,0)|U- 11/11 + |e| + |0| 

where ||.|| is defined by H2.4() and |.| is the absolute value in R. (|5.1() - (|5.2|I - H5.3|I - 
H5.4|l can be written in the following standard form of first order differential 
systems: 

\X{0) ^Xo \x{0)^{fo,eo,9o) 

We will prove the local existence theorem by applying the standard theorem for 
such systems in a Banach space, and that requires that F be continuous with 
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respect to t and locally Lipschitzian with respect to (/, e, 9) in the Banach space 
norm. In our case, the change of variable e = - shows that — y^l + a^jpP 
and the collision kernel A given by H3.4|l depends on e, So that, in (|5.1|I - H5.2|I - 
H5.3|l . Q, p, P depend only on f and e, and this implies that in H5.7|) . F does not 
depend explicitly on t, but only on f, e, 9 



5.2 The local Existence Theorem 



Proposition 5.1 There exist an interval [0,1], I > such that, the initial value 

problem for the system i|5. Jp - i|.5.ii)) - i|.5.5]) with initial data 

{fo, eo, ^o) G X R X M has a unique solution (/, e, 9) on [0, 1] 

Remark 5.1 For the proof of proposition 5.1, we will set 



Fi{eJ, 



1 



oQifJ), F,ieJ,9) 



= - — - 47r(p+3P)+A 
(5.8) 

in F2 and F^, (/, e, 0) 1-^ — |e and {f,e,9) 1-^ — ^ + A are C°° functions , which 
are then, locally Lipschitzian with respect to (f,e,9). The problem will be to 
prove that , in Fi and F3, (/, e) i-> j^Q(fJ) and (/, e) -47r(p + 3P)(/, e) 
are locally Lipschitzian with respect to (f e), in the L\{M!^) x M norm. 
We will deal with the quantities: 



1 



Q(/i,/i)(ei)- 



1 



P°(e2) 



Q(/2,/2)(e2); p(ei,/i)-p(e2,/2); P(ei, /i)-P(e2, /2) 



and this will lead us to look in detail inside the collision operator Q. It then 
shows useful to begin by giving some preliminary results. 

Lemma 5.1 We have for /, /i, /2 G Ll{M.'^), ei, 62 €]0, |], and C being a con- 
stant: 



1 



1 



P°(ei) pO(e2) 
Q(/i,/i)(e2)-g(/2,/2)(e2) 



Q{Lf){e2] 



< 



C 



l/ll|ei-e2| 



P°(e2) 

Proof of the lemma 5.1 

Recall that = 
1) we have: 

1 1 



<-5(ll/l||+ 11/211)11/2 -/ill 
^2 



(5.9) 
(5.10) 



p"(ei) P°(e2) 



But we have: 



1 + 
Q(/,,/)(e2) 



/(e2)-p''(ei) Q(/,/)(e2) 

X 



< 



pO(ei) 
1 



\P\' 



P°(e2) 



\P\' 



(p0(ei)+p0(e2))p0(ei) " {p^ie^W 



(P"(ei)+p0(e2))p0(ei) 
el^el<l 



Q(/,/)(e2) 



P°(e2) 



< 



\pI 2_ 2 
\P? 
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So we have, since < < 1 and using the expression 



of 



1 



1 



pO(ei) pO(e2) 



Q(/,/)(e2) 



< 



ere 



Q(/,/)(e2) 



P°(e2) 



|ei-e2| (a) 



l|5.9|l then follows from (a), the inequality H3.26|l on Q, in which we keep f, we 
set g = 0, and with a = ^ in H3.27|l . 

2) (|5.1U|) follows from (|3.2()|) in which we set: f — fi, g — f2 and a = — in 
llT?7ll . ■ 

Lemma 5.2 We have for f G LKR^), ei, 62 g]0, |] and C being a constant: 



[Q(/,/)(ei)-Q(/,/)(e2)] 



< 



ei - e2 



(5.11) 



Proof of Lemma 5.2 

Here we use Q = — , where Q^, Q~ are given by (|3.2I) and H3.3|l . in which 
a = i. We write: 

e 

/)(ei)-g(/, /)(e2) = [Q+(/, /)(ei)-Q+(/, /)(e2)] + [Q-(/, /)(e2)-Q-(/, /)(ei)] 

(a) 

We can, write, using the expression H3.2|l of 



[Q+(/,/)(ei)-Q+(/,/)(e2)] = 
where 



52 



i?(ei, e2,p, q')f{p')f{q')dqduo (b) 



B = B(ei,e2,p, q,p',q') 



1 



[A{e^yA{e2)V 



(c) 

in which A(ei) stands in fact for A{ei,p,q,p',q'), z = 1,2. Assumption 1)3. 7|l on 
the collision kernel A, gives: 



e?gO(ei) 



{A{e,) ~ Aie^)) 



< 



7 



e?<Z°(ei) 



ei - 62 



(d) 



Now we can write, using q^{ei) = . 1 + -^\q\'^ > and < Ci < 1, i — 1,2: 



e?gO(ei) eig"(e2) 



(ei-e?)gO(e2) + e?b"(e2)-<Z°(ei)] 



< 



9"(ei) 



(eie2)3g0(ei)g0(e2) 
3|ei - 62 



+ 



(6162)3 6?63(gO(62))2 



Hence: 



e?gO(6i) 63q0(e2) 



< 



elq^{ei) \e 



|6i - 62! 



(e) 
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(c) then gives, using (d), (e), assumption 1)3. 4|l on A, {\A\ < d), and < < 1: 
with C = 7 + 5Ci. We then deduce from (b) and (f), using expression 12 . 41 of 



,i [Q+(/i,/i)(ei)-Q+(/i,/i)(e2)] ,| ^ C|ei-e2| f f f VTTW\fm\fiQ')\ ,-,-, 

(g) 

We compute this integral the same way as in the proof of H3.24|l in lemma 3.1, 
using the change of variable {p, q) i— > (p', q') defined by (|3.11|) and this leads, 
using (|X77|) to 



[Q+(/i,/i)(ei)-Q+(/i,/i)(e2)] ^ C||/f|ei - 62! 



Next, we proceed the same way for the second term Q~{f, /)(e2) — /)(ei) 
of (a), using this time the expression (|3.3|) of Q~ . The only difference is that, in 
integrals (b) and (g) f{p')f{q') is replaced by f{p)f{q)- A direct computation 
without change of variable, then leads to the same estimation (h) where is 
replaced by Q", and lemma 5.2 follows. ■ 

Lemma 5.3 We have for ei,e2 S]0, |], /i,/2 G ^2(-'^'^) '^'^'^ C > being a 
constant. 



Q(/i,/i)(ei) Q(/2,/2)(e2) 



p"(ei) pO(e2) 

Proof of the lemma 5.3 

We can write: 



< g(ll/i|| + ll/iir + ll/2||)(|ei -e^l + Wh- /2I 



(5.12) 



Q(/i,/i)(ei) Q(/2,/2)(e2) Q(/i,/i)(ei)-Q(/i,/i)(e2) , / 1 1 



P°(ei) p0(e2) pO(ei) VP°(ei) P^ie^) 

+ ^jj^[0(/i,/i)(e2)-Q(/2,/2)(e2)] 

Then apply (|5.11|l with / = /i to the first term, I5.9|l with / = /i to the 
second term, and (|5.10() to the third term to obtain (|5.12|l . by addition of the 
inequalities and using < < 1: 

We now consider F3 in ()5.8|l . We can write since, by ()4.21(l and (|4.22|l . p and P 
are linear in f: 

P(ei, /i) - p{e2, /2) = p(ei, /i - /2) + [p(ei, /2) ~ p(e2, /2)] (5.13) 
P(ei, /i) - P(e2, /2) - P(ei, /i - /2) + [P(ei, /2) - P(e2, /2)] (5.14) 



Q(/i,/i)(e2) 
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gives, using Jl + m^< \^T+W 



|p(ei,/i-/2)| < ^ll/i-MI (5.15) 
H5.15|l gives, using (|4.23() (i.e P < p): 

|P(ei,/i-/2)| < ^||/i-/2|| (5.16) 

We then deduce from (|5.13() and H5.14|l . using (|4.37fl . (|4.38|l in which we set 
/ = /2, (EEl), and < ei < 1, i = 1 , 2, that 

|p(ei. A) - p(e2, /2)| < ^^^^^^(|e2 - ei| + ||/i - Ml) (5.17) 

^1^2 

|P(ei,/i)-P(e2,/2)| <^^^^|^(|e2-ei| + ||/i-/2||) (5.18) 

Now let /I G L^(M^) and g]0, §]; take ei,e2 and 

A, /2 G B{f\ 1) := {/ G L\{m.^)\\f - /I < 1|} 

Then we have , J-, i < |- and ||/i||, UMI < + 1. 

We then deduce from ifCTTjl and that: 



Q(/i,/i)(ei) Q(/2,/2)(e2) 



<^^(ll/i-/2|| + |ei-e2|) 



P°(ei) p"(e2) 

|p(ei,/i) - p(e2, /2)| < K(\\h - h\\ + |ei - 62!) 

|P(ei,/i) - P(e2, h)\ < K{\\h - /2II + |ei - 62!) 

where K — K{e^, f^) is a constant depending only on and If we add this 
to the fact that {f,e,6) 1-^ — |e and {f,e,6) 1-^ — ^ + A are Lipschitzian with 
respect to the norm of P = L2(M'^)xlRxR, we can conclude that F = (Pi, P2, P3) 
defined by (|5.9|l is locally Lipschitzian in {f,e,9) with respect to the norm 
of E. Proposition 5.1 then follows from the standard theorem on first order 
differential system for functions with values in a Banach space. Notice the e, 9 
are continuous and / G C[0, /; Ll{R^)] ■ 
From proposition 5.1, we deduce the following theorem: 

Theorem 5.1 Let /o G i2(^^); /o > a.e, ap > § and let a strictly positive 
cosmological constant A, be given. Define ciq by the relation ^5.6]) . 
Let r > |j/o|| . Then, there exist a number I > such that the initial value 
problem for the coupled Einstein- Boltzmann system \2.1!^) - \4--l^ - {4--^ with the 
the initial data {fo, 00,0-0), has an unique solution (f a) on [0,1]. The solution 
(a, f) has the following properties: 

(i) a is an increasing function. 
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(ii) 

(in) 



f eC[[Q,l],Xr] (5.19) 
Ill/Ill <ll/o|| (5.20) 



Proof of theorem 5.1 

Choose in proposition 5.1, /o,eo and as in (|5.5|) . Let {e,f,9) be the unique 
sohition of that initial value problem defined on the interval [0, /], / > 0, whose 
existence is proved by that proposition. Set a — ^, 1)5.2(1 then implies 6 = 3^. 
We know, by the characteristic method that the I3oltzmann equation 12.13|l is 
equivalent to 1)5. 1|) . We also know that, in the framework described in §4.4, and 
in which we solve the Einstein equation in the case A > 0, the system (14. If) - (14.211 
in a, and the system H5.2|l - H5.3|l in (e, 9) are equivalent. Since the relation (|5.6|l 
which implies (|4.1Q(I . is satisfied, we know that the Hamiltonian constraint (|4.1)l 
is satisfied on [0,^]. We can conclude that, in the framework fixed by choosing 
1)4.16)1 in the case A > , the unique solution {f,e,9) on [0,^], of the initial 
value problem ()5.1)) - (|F?^ - I)5.3)I - ()5.4)I with /o, eg, 6'ogiven by 1)5.5)) and satisfying 
1)5. 6)) . give the unique solution (/, a — i) on [0,1], of the initial value problem 
I)2.13)) - I)4.1)) - (I4.2)) with the initial data (/o,ao,ao) satisfying ()5.6)) . 
Concerning the properties of the solution (f, a) on [0,^[: ()4.16)) shows that 
a > 0; i I— > a{t) then satisfies the point i) and a is bounded from below, since 
a{t) > ao > |. a then satisfies the hypotheses of theorem 3.1 in which we set 
to^O,T = 1. Note that C[[0,l];Xr] C C[[0,l]; Ll{R^)]; fTT^ follows from the 
uniqueness and 1)5.20)) from 1)3. 18)) . B 

6 Global Existence Theorem for the Coupled 
Einstein-Boltzmann System 

6.1 The Method 

Here, we prove that the local solution obtained in §5 is, in fact, a global solution. 
Let us sketch the method we adopt: Denote [0,T[, where T > 0, the maximal 
existence domain of the solution of ()5.1)) - ()5.2)) - l)5.3)l . with initial data (/o, eg, ^o) 
defined by ()5.4)) - ()5.5)) - ()5.t))l : here we denote this solution (/, e, 9), in order words 
we have, on [0, T[ 

/ = ^W-/) (6.1) 
e = --e (6.2) 

n2 

9 = -—-47:{p + 3P)+A (6.3) 

O 

/(0) = /oeX,, g(0)^eo--, m^3^ (6.4) 

ao ao 
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with in Ht).3|l p = p{f,e): P = P{f,e) and in Ht).4|l /o, aq, 6*0 subject to the 
constraints (|5.6() . r > is given such that r > ||/o||. 

If T = +00, the problem is solved. We are going to show that, if we suppose 
that T < +00, then the solution {f,e,6) can be extended beyond T, which 
contradicts the maximality of T. Suppose < T < +00 and let to € [0, T[. We 
will show that, there exists a strictly positive number S > 0, independent of 
to, such that the following system in (e, 6) on [to, to + 6], in which a 



1 . 



/ = (6.5) 



e = --e (6.6) 

n2 

^ = -— -47r(p + 3P)+A (6.7) 
o 

/(io)=/>o), e(to) = e~(to), 9 (to) = 9 (to) = 3^ (6.8) 

a (to) 

has a solution {f,e,6) on [to, to + 5]- Then, by taking to sufhciently close to T, 
for example, to such that < T — to < |, hence T < to + |, we can extend 
the solution {f,e,9) to [0, to + |[ that contains strictly [0,T[, and this will 
contradicts the maximality of T.We need some preliminaries results. 
In what follows, we suppose < T < +00 and to G [0, T[. 

6.2 The Functional Framework 

In all what follows, C2, C3 , Dq are the absolute constant defined by (|4.36|) and 
(II3SJ1. We set, for 6 > 0: 

< - |eeqto,to + ^],^e-^^(*«+*+i)' <e(to+t)< ^ ,Vte [0,^[| 



{0 G C[to, to + S],0, V3A < 6i(to + t) < Do Vt e [0, 6[} 



where C[to, to + S] is the space of continuous (and hence bounded) functions on 

[to, to + S]. One verifies easily that Ef^^ and Fj"^ are complete metric subspaces 

of the Banach space (C[to, to + 5], ||.||oo) where ||m||oo = Sup \u{t)\ 

te[to4o+5] 

6.3 The Global Existence Theorem 

Proposition 6.1 There exists a strictly positive real number S > depending 
only on the absolute constants oq, A , r and T such that the initial value problem 
has a solution (/, e ^,6) G C[[to,to + S]; Xr]x Ef^x F^^ 
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Proof 

It will be enough, if we look for 5 such that < 5 < 1. By theorem 3.1, we know 
that if we fix e S E^^^ and if we set a = 4, then IjCSfl has an unique solution 
/ e C[[to,tn + 5]\Xr\, such that, /(to) = /(to), and, by (jXl^ and 



Wfm < ||/>o)|| < ll/oll <r 



(6.9) 



Next, by proposition 4.6 in which we set :z — e,Q — 9, we know that if / 
is given in C[[to!^o + <5];-'^r]i then H6.6|) - (|6.7|) has an unique solution {e,9) on 
[to, to + 6] such that e(to) = 6'(to) = 3|||^. Now and (|07l) show 

that (e ~ ■ 



e Ef^ X Ff^_^ . This allows us to define the application 



G ■.C[[to,to + 5];Xr]x El 



C[[to,to + 5];Xr] X « xF/j (6.10) 
G(/,e-) = [/,(e,0)] (6.11) 



We are going to show that we can find (5 > such that G defined by 16.1()|l 
induces a contracting map of the complete metric space C[[to, to + S];Xr] x Ef^ 
into itself, that will hence, have an unique fixed point (/, e) ; this will allow 
us to find 9 such that (/, e, 9) be the unique solution of 



in C[[to,to + 6];Xr] x (£;*, x F°J. So if we set in 



^ e = e e E, 



to , 



P = P = pie J), P = P = Pie J) where / G C[[to,to + S];Xr], we have 
a solution (/, e, 9) of that system, or, equivalently, a solution (/, e, 0) of the 
following integral system: 



fito + t) = /(to) + 



to+t 



1 



to 



pO(e-) 



QifJ)ie)is)ds 



e(to + t) = e(to) + 



to-t 



9is)e{s] 



ds 



9{to + t)^9{to) + 



to+t 



to 



[---4^(p + 3P) 



A](s)ds 



(6.12) 

(6.13) 
(6.14) 



tG[0,,5[ 

Let 61,62 e El^;fi,f2 e C[[to, to+(5]; -Xr]. To (r esp /» ) i ^ 1 , 2, corresponds 
by G, the solution /^ (resp (6^,6'^)) of l|0^ . [resp (|F!ni) - (|^?TII) ]. Writing each 
equation for i = 1, i = 2 and subtracting yields: 



ll/l-/2||< 



to+t 



to 



pO(ei) 



Q(/i,/i)(ei)- 



1 



P°(e2) 



Q(/2,/2)(e2) 



jei - 62! < 



1^1-^2! < 



to+t 



to 



ds (6.15) 
(6.16) 

9'i\ + 127r(|pi - /J2I + [A - P2\)]ds (6.17) 



rto+t 


9iei 


^'262 


'to 


3 


3 



ds 
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We are going to bound the r.h.s of Ht).15|l and (|6.17ll . using: l|5.12|) in which 
we set ei = ei, = e.2 and (|5.17|I - H5.18|I . in which we set /i — /i, /2 = /2. 
Since 61,62,61,62 € i^j^g, (jO^ . shows that: < €26^''^^+'^'^^ , since to < T, 

^ < S < 1_. So we deduce from (6.12), (6.13)' and (6.14), using (jO^ . (|CT7|) . 
(ElSJI/dll <r ,0 <t < S, th e definition of i;^*^, (|6'i| < Do{ao,r, A,T)), 
see ()4.45|) . the definition (|2.16|l of the norm |||.||| of function over [to, to + 6], 
that: 

III/1-/2III <^Mi(|lei-e2||oo + |||/i-/2|||) (6.18) 



|ei - 62II00 < 6M2(\\ei - 62! 



dloo) 



(6.19) 



11^1 - 021100 < SAhiWei - 62IU + 11^1 - 02\\oo + II 1/1 - /2III 



(6.20) 



where Mi, M2, M3 are constants depending only on ao. A, r and T. We have by 
addition of KT^ and (|OII) : 

||6l-62||oo + ||^?l-^?2||oc.) <2(5(Af2+M3)(|l6i-e2||oo + ||01-02||oo + |||,/l-/2|||) 

(6.21) 

Then, if we choose S such that: 



S = Inf 



1 



8(Mi +M2 +M3)_ 
Ij^:^ imphes that : SAh < i; 26{M2 + Afg) < i and ^J^, (??) give: 

lll/i - /2III < ^l|ei - e2||oo 



(6.22) 



[||6l-e2|U + 11^1 -021100 <i|ll/l-/2|l 



and by addition: 



lll/i - /2III + Il6i - 62II00 + 1101 - 02IIOO < 3(111/1 - Ml! + llei - e2|loo) 
from which, we deduce: 



fl-/2|l| + |l6l-62|loo < 3(111/1-/2111 + 1161-621100) 



(6.23) 



(/, e) is a contracting map from the complete 



H6.23|l shows that the map (/, e 
metric space C[[to , to +^] ; Xr] x Ef^^ into itself, for every S satisfying (|6.22|) , which 
shows that such a 6 depends only on ao,r, A and T. This map has an unique 
fixed point (f, e); since e is known, H6.6|l determines hy 9 ^ —3-, and (/, 6, 6) 
is a solution of (ESJ-ESl-EiZIl-EiHl) in C[[to,to + S];Xr] x £;£ x i^/^. This 
completes the proof of proposition 6.4 ■ 
We can then state: 
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Theorem 6.1 The initial value problem for the Einstein- Boltzmann system 
with a strictly positive cosmological constant K on a Robertson- Walker space- 
time has a global solution (a, f) on [0, +oo[, for arbitrarily large initial data oq 
and fo e Ll{R^), /o > a.e . 

Remark 6.1 1) Nowhere in the proof we had to restrict the size of the initial 
data Oq, fo, which can then be taken arbitrarily large. 

2) In the author considered only the Hamiltonian constraint ^4-1^ , in 
the case A = 0, without studying the evolution equation i4.'J^ that is, as 
we saw, the main problem to solve, since the Hamiltonian constraint 

is satisfied once it is the case for the initial data. 

3) We will prove in a future paper, that theorem 6.1 extends to the case A = 

4) In the future, we will try to relax hypotheses on the collision kernel A. 
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